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Abstract 

The production of Z° bosons in the reaction ep — > eZ°p(*\ where p^ stands 
for a proton or a low-mass nucleon resonance, has been studied in ep collisions 
at HERA using the ZEUS detector. The analysis is based on a data sample 
collected between 1996 and 2007, amounting to 496 pb -1 of integrated lumin- 
osity. The Z° was measured in the hadronic decay mode. The elasticity of 
the events was ensured by a cut on ?7 max < 3.0, where ?7 max is the maximum 
pseudorapidity of energy deposits in the calorimeter defined with respect to 
the proton beam direction. A signal was observed at the Z° mass. The cross 
section of the reaction ep — > eZ°p(*> was measured to be a [ep — > eZ°p ( -*^ = 
0.13 ± 0.06 (stat.) ± 0.01 (syst.) pb, in agreement with the Standard Model pre- 
diction of 0.16 pb. This is the first measurement of Z° production in ep collisions. 
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1 Introduction 



The production of electroweak bosons in ep collisions is a good benchmark process for test- 
ing the Standard Model (SM). Even though the expected numbers of events for W ± and 
Z° production are low, the measurement of the cross sections of these processes is import- 
ant as some extensions of the SM predict anomalous couplings and thus changes in these 
cross sections. A measurement of the cross section for W ± production at HERA has been 
performed by HI and ZEUS pQ in events containing an isolated lepton and missing trans- 
verse momentum, giving a cross section a {ep — » W ± X) = 1.06 ± 0.17 (stat. © syst.) pb, in 
good agreement with the SM prediction. The cross section for Z° production is predicted 
to be 0.4 pb. 

This paper reports on a measurement of the production of Z° bosons in e ± p collisions 
using an integrated luminosity of about 0.5 fb -1 . The hadronic decay mode was chosen^ 
because of its large branching ratio and because it allows the excellent resolution of the 
ZEUS hadronic calorimeter to be exploited to the full. The analysis was restricted to 
elastic and quasi-elastic Z Q production in order to suppress QCD multi-jet background. 
The selected process is ep — > eZ°p(*\ where p^ stands for a proton (elastic process) or a 
low-mass nucleon resonance (quasi-elastic process). 

Figure [T] shows a leading-order (LO) diagram of Z° production with subsequent hadronic 
decay. In such events, there are at least two hadronic jets with high transverse energies, 
and no hadronic energy deposits around the forward^! direction, in contrast to what would 
be expected in inelastic collisions. 



2 Experimental set-up 

HERA was the world's only high-energy ep collider, with an electron^] beam of 27.6 GeV 
and a proton beam of 920 GeV (820 GeV until 1997). For this analysis, e ± p collision 
data collected with the ZEUS detector between 1996 and 2007, amounting to 496 pb -1 of 
integrated luminosity, have been used. They consist of 289 pb -1 of e + p data and 207 pb -1 
of e~p data. 

1 The Z° decay into charged lepton pairs was studied in a previous ZEUS publication [2]. 

2 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the forward direction, and the X axis pointing towards the 
centre of HERA. The coordinate origin is at the nominal interaction point. The pseudorapidity is 
defined as rj = — In (tan |), where the polar angle, 9, is measured with respect to the proton beam 
direction. 

3 The term electron also refers to positrons if not stated otherwise. 
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After 2003, HERA was operated with a polarised lepton beam. When combining the data 
taken with negative and positive polarisations, the average polarisation is less than 1% 
and its effect was neglected in this analysis. 

A detailed description of the ZEUS detector can be found elsewhere [3]. A brief outline 
of the components that are most relevant for this analysis is given below. 

Charged particles were tracked in the central tracking detector (CTD) [4], which operated 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The CTD 
consisted of 72 cylindrical drift chamber layers, organised in nine superlayers covering 
the polar-angle region 15° < 9 < 164°. For the data taken after 2001, the CTD was 
complemented by a silicon microvertex detector (MVD) [5j, consisting of three active 
layers in the barrel and four disks in the forward region. 

The high-resolution uranium-scintillator calorimeter (CAL) [6] consisted of three parts: 
the forward (FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part was 
subdivided transversely into towers and longitudinally into one electromagnetic section 
(EMC) and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). 
The smallest subdivision of the calorimeter was called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, were a(E)/E = 0.18/ y/E for electrons and 
cr(E)/E = 0.35/ VE for hadrons, with E in GeV. 

The luminosity was measured using the Bethe-Heitler reaction ep — > cjp by a luminos- 
ity detector which consisted of independent lead-scintillator calorimeter [7] and magnetic 
spectrometer [8] systems. The fractional systematic uncertainty on the measured lumin- 
osity was 2 %. 

3 Monte Carlo simulations 

Monte Carlo (MC) simulations were made to simulate the Z° production process. They 
were used to correct for instrumental effects and selection acceptance and to provide a 
template for the shape of the invariant- mass distribution of the Z° signal. The EPVEC 
program [9] was used to generate the signal events at the parton level. The following Z° 
production processes are considered in EPVEC: 

• elastic scattering, ep — > eZ°p, where the proton stays intact; 

• quasi-elastic scattering, ep — > eZ°p*, where the proton is transformed into a nucleon 
resonance p*; 

• deep inelastic scattering (DIS), 7*p — > Z°X, in the region Q 2 > 4 GeV 2 , where Q 2 is 
the virtuality of the photon exchanged between the electron and proton; 
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• resolved photoproduction, — > (qq — > Z°) X, where one of the quarks is a constituent 
of the resolved photon and the other quark is a constituent of the proton. 

In EPVEC the first two processes are calculated using form factors and structure functions 
fitted directly to experimental data. Note that, even if the virtuality of the exchanged 
photon is small, the scattered electron could receive a large momentum transfer when 
the Z° is radiated from the lepton line. In the last two processes, the proton breaks up. 
The DIS process is calculated in the quark-parton model using a full set of leading-order 
Feynman diagrams. Resolved photoproduction is parametrised using a photon structure 
function and is carefully matched to the DIS region. The cross section of Z° production 
is calculated to be 0.16 pb for elastic and quasi-elastic processes and 0.24 pb for DIS 
and resolved photoproduction. The difference between e + p and e~p cross sections is 
negligible for this analysis (jl% for the DIS process). A correction, based on the MC 
cross section, was made to account for the part of data taken at the centre-of-mass energy 
y/s = 300 GeV, so that the result is quoted at y/s = 318 GeV. 

After the parton- level generation by EPVEC, PYTHIA 5.6 [TU] was used to simulate 
initial- and final-state parton showers with the fragmentation into hadrons using the Lund 
string model [11] as implemented in JETSET 7.3 [10]. The generated MC events were 
passed through the ZEUS detector and trigger simulation programs based on GEANT 
3.13 |12] . They were reconstructed and analysed by the same programs as the data. 

A reliable prediction of background events with the signal topology, which are predom- 
inantly due to the diffractive photoproduction of jets of high transverse momentum, is 
currently not available. Therefore, the background shape of the invariant-mass distribu- 
tion was estimated with a data-driven method. The normalisation was determined by a 
fit to the data. 

4 Event reconstruction and selection 

The events used in this analysis were selected online by the ZEUS three-level trigger 
system [13], using a combination of several trigger chains which were mainly based on re- 
quirements of large transverse energy deposited in the calorimeter. In the offline selection, 
further criteria were imposed in order to separate the signal from the background. 

The events are characterised by the presence of at least two jets of high transverse energy 
and, for a fraction of events, by the presence of a reconstructed scattered electron. In order 
to select events with a Z° decaying hadronically, jets were reconstructed in the hadronic 
final state using the kx cluster algorithm [TJ] in the longitudinally invariant inclusive 
mode [15]. The algorithm was applied to the energy clusters in the CAL after excluding 
those associated with the scattered-electron candidate [16TfTB"]. Energy corrections [T§H2"T] 
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were applied to the jets in order to compensate for energy losses in the inactive material 
in front of the CAL. 

In this analysis, only jets with Et > 4GeV and \r)\ < 2.0 were used. Here Et is the jet 
transverse energy and r\ its pseudorapidity. The hadronic Z° decay sample was selected 
by the following requirements on the reconstructed jets: 

• at least two jets in the event had to satisfy Et > 25 GeV; 

• |A0j| > 2rad, where A0j is the azimuthal difference between the first and second 
highest-E'T jet, as the two leading jets from the Z° boson decays are expected to be 
nearly back-to-back in the X—Y plane. 

Electrons were reconstructed using an algorithm that combined information from clusters 
of energy deposits in the CAL and from tracks [TB]. To be defined as well-reconstructed 
electrons, the candidates were required to satisfy the following selection: 

• E' e > 5 GeV and E in < 3 GeV, where E e is the scattered electron energy and E m is the 
total energy in all CAL cells not associated with the cluster of the electron but lying 
within a cone in r] and of radius R = a/ A?] 2 + A0 2 = 0.8, centred on the cluster; 

• If the electron was in the acceptance region of the tracking system, a matched track 
was required with momentum ptrack > 3 GeV. After extrapolating the track to the 
CAL surface, its distance of closest approach (DCA) to the electron cluster had to be 
within 10 cm. 

The following cuts were applied to suppress low-Q 2 neutral-current and direct-photopro- 
duction backgrounds: 

• -Ercal < 2 GeV, where -Ercal is the total energy deposit in RCAL; 

• 50 < E — pz < 64 GeV, where E —pz = ^ Ei (1 — cos 0$); Ei is the energy of the z-th 
CAL cell, 6i is its polar angle and the sum runs over all cellsj; 

• 9 e < 80° for well reconstructed electrons, where 9 e is the polar angle of the scattered 
electron, motivated by the fact that, due to the large mass of the produced system, 
the electron is backscattered to the forward calorimeter or forward beam pipe; 

• the event was rejected if more than one electron candidate was found; 

• jets were regarded as a misidentified electron or photon and were discarded from the 
list of jets if the direction of the jet candidate was matched within R < 1.0 with that 
of an electron candidate identified by looser criteria^ than those described above. This 

4 For fully contained events, or events in which the particles escape only in the forward beam pipe, the 
E — pz value peaks around twice the electron beam energy, 55 GeV. 

5 Candidates were selected by less stringent requirements and clusters with no tracks were also accepted 
to find photons and electrons. 
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cut causes a loss of acceptance of about 3%. 

To remove cosmic and beam-gas backgrounds, events fulfilling any of the conditions listed 
below were rejected: 

• \Z vtx \ > 50 cm, where Z vtx is the Z position of the primary vertex reconstructed from 
CTD+MVD tracks; 

• 175° < (#jeti + 0jet2) < 185° and Acfij > 175° simultaneously, where 9j etl and 6- ]et 2 are 
the polar angles of the first and second highest- .Et jet, respectively, and A<ftj is the 
azimuthal difference between them; 

• \t u — td\ > 6.0 ns, where \t u — t d \ is the timing difference between the upper and the 
lower halves of the BCAL; 

• j>T > 25 GeV, where j>T is the missing transverse momentum calculated from the 
energy clusters in the CAL; 

• < 0.25 (iV^ — 20), where jVJjjf is the number of tracks associated with the 
primary vertex and is the total number of tracks [22]- 

The number of events passing the above selection was 5257. Finally, to select the elastic 
and quasi-elastic processes preferentially, a cut on r] max was introduced, 

• ??max < 3.0. 

The quantity ^ max was defined as the pseudorapidity of the energy deposit in the calori- 
meter closest to the proton beam direction with energy greater than 400 MeV as determ- 
ined by calorimeter cells. This cut also rejected signal events which have energy deposits 
from the scattered electron in the calorimeter around the forward beam pipe, causing an 
acceptance loss of about 30%. 

After all selection cuts, 54 events remained. The total selection efficiency was estimated 
by the MC simulation to be 22% for elastic and quasi-elastic processes and less than 1% 
for DIS and resolved photoproduction events. The number of expected signal events in 
the final sample, as predicted by EPVEC, is 18.3. The contribution from elastic and 
quasi-elastic processes amounts to 17.9 events. 

5 Background-shape study 

Figure [2^ shows the distribution of the invariant mass, Mj et s, after all the selection cri- 
teria except for the requirement 77 max < 3.0. The variable Mj ets was calculated using all 
jets passing the selection criteria described in SectionHJ Figures |2b-d show Mj ets for vari- 
ous r] max slices in the inelastic region (r7 max > 3.0) for the same selection. No significant 
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dependence on r] mSLX of the Mj e t s distribution beyond that expected from statistical fluctu- 
ations was observed in the inelastic region. In addition, the shape of the Mj ets distribution 
outside the Z° mass window in the region r/ max < 3.0 was found to be consistent with that 
in the inelastic region (Fig. [3]). Therefore, the Mj ets distribution in the inelastic region was 
adopted as a background template in a fit to the data in the elastic region as described 
in the following section. 



6 Cross-section extraction 

A fit to the sum of the signal and a background template for the Mj ets distribution was 
used for the cross-section extraction. The template A^y is defined according to: 

N rei , = aN$l (e) + bN^f, (1) 

where i is the bin number of the Mj ets distribution. The parameter e accounts for a pos- 
sible energy shift, i.e. Mj ets = (1 + e) M^, where Mj^ is the invariant-mass distribution 
of the signal Z° MC. The quantity JVjJJ 3 is a signal template estimated from the Z° MC 
distribution after all cuts, normalised to data luminosity. The quantity A^|* a is a back- 
ground template determined from the data outside the selected region. The parameters 
a and b are the normalisation factors for the signal and background, respectively. The 
likelihood of the fit, C, is defined as follows: 

C = C 1 (N ohs ,N mi )xC 2 {e,a e ), (2) 

with 

T-r exp(-iV refi ) (A/ re fi) JVobM , „ ft 2 

Here £i (N ^ s , A^ re f) is the product of Poisson probabilities to observe AT bs,j events for 
the bin i when N Te fj is expected. The term £ 2 (e, &e) represents the Gaussian probability 
density for a shift e of the jet energy scale from the nominal scale, which has a known 
systematic uncertainty of a e = 3%. From the likelihood, a chi-squared function is defined 

as 



l~\ l^Vob s ,^Vobs) ^ \cr e J 



with 



N ieiji - N ohs>i + N ohs , t ln (N ohSji /N iei>i ) (if N ohStl > 0) 

N le{>i (if N ohSii = 0) . 

The best combination of (a,b,e) is found by minimising x 2 ■ The value of a after this 
optimisation gives the ratio between the observed and expected cross section, i.e. <7 b s = 
ao"sM- The maximum and minimum values of a in the interval Ax 2 < 1 define the range 
of statistical uncertainty. 
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7 Systematic uncertainties 



Several sources of systematic uncertainties were considered and their impact on the meas- 
urement estimated. 

• An uncertainty of 3% was assigned to the energy scale of the jets and the effect on 
the acceptance correction was estimated using the signal MC. The uncertainty on the 
Z° cross-section measurement was estimated to be +2.1% and —1.7%. 

• The uncertainty associated with the elastic and quasi-elastic selection was considered. 
In a control sample of diffractive DIS candidate events, the ?? max distribution of the 
MC agreed with the data to within a shift of ?7 max of 0.2 units [23]. Thus, the r/ max 
threshold was changed in the signal MC by ±0.2, and variations of the acceptance 
were calculated accordingly. The uncertainty on the cross-section measurement was 
+6.4% and -5.4%. 

• The background shape uncertainty was estimated by using different slices of ?7 max in 
the fit. The background shape was obtained using only the regions of 4.0 < r/ max < 4.2 
or 4.2 < i] max . The region of 3.0 < i] max < 4.0 was not used since a small number 
of signal events is expected in this r] max region§. The resulting uncertainty in the 
cross-section measurement was ±1.5%. 

• The uncertainty associated with the luminosity measurement was estimated to be 2%, 
as described in Section [2J 

• The Z° mass distribution from the MC used as a signal template has a Gaussian core 
width of 6 GeV. A possible systematic uncertainty coming from the width of the 
MC signal peak was studied. The mass fit was repeated after smearing the Z° mass 
distribution in the MC by a Gaussian function with different widths. The measured 
cross section did not change significantly after smearing the distribution up to the 
point where the fit x 2 changed by 1. No systematic uncertainty from this source was 
assigned. 

The total systematic uncertainty was calculated by summing the individual uncertainties 
in quadrature and amounts to +7.2% and —6.2%. 

8 Results and conclusions 

Figure [3] shows the invariant-mass distribution of the selected events. It also shows the 
fit result for the signal plus background and the background separately. The fit yielded a 

6 The ratio of the expected number of signal MC events to the observed data in this region was estimated 
to be 2.6% for 80 < Afj Cts < 100 GeV, while in the other slices it was less than 0.4%. 
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result for the parameter a from Eq. [TJof a = 0. 821q35- That translates into a number of 
observed Z° events of 15.0^;° (stat.), which corresponds to a signal with a 2.3 a statistical 
significance. The fit yielded a value for the energy shift e of 0.028^q 020 > which is compatible 
with zero. The quality was evaluated according to Eq. 3; the value of x 2 /ndf = 17.6/22, 
where ndf is the number of degrees of freedom, indicates a good fit. The cross section for 
the elastic and quasi-elastic production of Z° bosons, ep — > eZ°p(*\ at yfs = 318 GeV, 
was calculated to be 

a(ep -+ eZ°p M ) = 0.13 ± 0.06 (stat.) ± 0.01 (syst.) pb. (4) 

This result is consistent with the SM cross section calculated with EPVEC of 0.16 pb. 
This represents the first observation of Z° production in ep collisions. 
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Figure 2: The Mjets distribution of the data (a) after all selection criteria, except 
for the rj msix cut, (b-d) in several rj max slices. 
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Figure 3: The Mj ets distribution and the fit result. The data are shown as points, 
and the fitting result of signal+background (background component) is shown as 
solid (dashed) line. The signal contribution is also indicated by the shaded area and 
amounts to a total number of N t, s events. The error bars represent the approximate 
Poissonian 68% CL intervals, calculated as ±y/n + 0.25 + 0.5 for a given entry n. 
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